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We present results for the QCD Equation of State (EoS) obtained using simulations of lattice 
QCD at zero chemical potential. Our high statistics results compare improved asqtad and p4fat3 
staggered quarks on lattices with a temporal extent N z = 6 and 8 and light quark masses approxi- 
mately one fifth and one tenth the strange quark mass. We find that the two actions give consistent 
results and estimate that the trace anomaly (e — 3p)/T 4 obtained on N T = 8 lattices represents the 
continuum value to better than 20% uncertainty over the temperature range 140 — 700 MeV. The 
precision in the estimates of energy density and pressure are better, therefore, we conclude that 
lattice estimates of the energy density and pressure should be used in the phenomenological anal- 
ysis of RHIC and LHC data. We also find a consistent picture of the crossover temperature from 
all observables studied, with the best estimated range to be 185 — 195 MeV. These calculations 
are being carried out on the IBM BlueGene/L supercomputer at Lawrence Livermore National 
Laboratory and at the New York Center for Computational Science (NYBlue). 
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1. Introduction 

Experiments at the Brookhaven Relativistic Heavy Ion Collider (RHIC), and at the LHC, will 
study the formation and evolution of the quark gluon plasma. The phenomenological analysis 
of the data has relied on a hydrodynamic description of the medium, and the results suggest a 
strongly interacting plasma with very small viscosity. A crucial input (and eventually the goal) 
of this phenomenological analysis is the characterization of the equation of state (EoS) over the 
temperature range 140 — 700 MeV relevant to experiments at both RHIC and LHC. Theoretically, 
simulations of Lattice QCD provide a first principles analysis of QCD in the vicinity of thermal 
equilibrium and zero baryon and strangeness chemical potential (For recent reviews see [|l],[2|]). The 
results of our detailed analysis of the equation of state as a function of temperature will, therefore, 
provide crucial guidance in the phenomenological interpretation of experimental measurements. 

A related quantity of high interest is the transition temperature from hadrons to a quark-gluon 
plasma (QGP). The minimum energy density required to produce a quark-gluon plasma grows as 
the fourth power of the temperature. Thus, a 10% error in the threshold temperature corresponds 
to a 45% error in the threshold energy density. Previous calculations with staggered fermions have 
found the transition for 2 + 1 flavors to be a rapid crossover. In light of this lack of a true phase 
transition (see, for example, []3|]), we will discuss what quantities to focus on for estimating the 
transition temperature needed in the phenomenological analysis of heavy ion collisions. 

This talk updates the status of HotQCD results for the EoS and the transition temperature ob- 
tained from large scale simulations on the BlueGene L at LLNL and NYCCS (NYBlue) [|, 0, g]. 
These calculations have been carried out with two sets of €?(a 2 ) improved actions, asqtad and 
p4fat3 staggered fermions, and for each of these actions we have simulated on lattices with extent 
Nt = 6 and 8 in the Euclidian time direction and a spatial extent of 32 3 . The strange quark mass is 
fixed to roughly the physical strange quark value, and mi/m s = 0.1 and 0.2 corresponding to Gold- 
stone M n 215 and 304 MeV. For preliminary results at m,( = 0.05m v see talk by Soeldner [§]. The 
associated T = calculations needed to subtract ultraviolet divergences in the EoS and determine 
the lattice spacing were done on 32 4 or larger lattices. We discuss uncertainties associated with the 
continuum and chiral extrapolations based on these combined data sets. 

2. Parameter Sets Used in the Simulations 

The run parameters for simulations with the p4fat3 action are given in Table [IJ, and those for the 
asqtad action in Table ||. The lattice scale is set using r$ (or equivalently r\ as they give consistent 
estimates). The simulations are being carried out along lines of constant physics by adjusting the 
bare strange quark mass to produce an approximately constant physical value of M Ss = 686 MeV 
corresponding to M^ro = L58 and M n r§ fh 0.52 for Me = 0.lm s along the trajectory. Having fixed 
the strange quark mass we are simulating three light quark mass values held fixed at mi/m s = 
0.2, 0.1 and 0.05. To allow comparison with previous studies the asqtad trajectory in HotQCD 
simulations is set 20% higher than the physical strange quark mass. 

3. Strategy for Precision Calculations 

Our goal is to map out the EoS over the temperature range 140-700 MeV to within 5% uncer- 
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Table 1: Simulation parameters for the p4fat3 action. For thermalization 800 trajectories are discarded. We 
also indicate the number of independent streams used to accumulate statistics. 
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Table 2: Simulation parameters for the asqtad action. 1000-1200 time units are discarded for thermalization. 

tainty. In addition to performing high statistics simulations (see Tables |j] and ||) we have adopted 
the following strategy to understand and control systematic errors. 

1. Discretization Errors: We are simulating two improved staggered actions - asqtad and p4fat3 
- that have different 0(a 2 ) errors. This provides a check but does not address the rooting 
issue for staggered fermions. 

2. Continuum Limit: Our goal is to perform a continuum extrapolation along lines of constant 
physics using lattices with N t = 6, 8 and 12. 

3. Extrapolation to physical u,d quark masses and the chiral limit will be done using mtjm s = 
0.2, 0.1 and 0.05. Staggered taste violations in finite T calculations require resolving what 
effective M n should correspond to these quark masses and thus the physical value of m u j [jj]]. 

4. Crossover Temperature: Recognizing the absence of a phase transition at physical values of 
the quark masses, we are simulating at 2 — 5 MeV interval over the entire crossover region 
to provide a precise quantitative picture of the transition in energy density, pressure, etc. 

5. T = simulations needed for performing subtractions of lattice artifacts in the determination 
of the EoS: We are simulating almost as many j3 values (« 20 for each action at N t = 8) 
as used in finite temperature runs. Estimates of a includes the much larger asqtad zero- 
temperature program. 
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4. Equation of State 

Our preliminary results for the trace anomaly ((e — 3p)/T 4 ) for both actions and for N t = 6 
and 8 lattices are shown in Fig. [j] Over the full range T = 140 — 500 MeV, we find the same pattern 
for both actions on going from N t = 6 to 8, i.e. a decrease in peak height and a slight shift of points 
to lower T, with the largest change in the range T = 190 — 300 MeV. A more detailed picture of 
the data over (140 < T < 200) MeV is also presented in Fig. [| (right figure), and in Fig. ^ for 
(180 < T < 300) MeV (left figure), and (300 < T < 700) MeV (right figure). 
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Figure 1: Preliminary results for the trace anomaly for the p4fat3 and asqtad actions as a function of tem- 
perature in MeV and in units of r$. The figure on the right magnifies the range T = 140 — 200 MeV and 
includes fits to the p4fat3 points. The purple curve above the data is from the hadron resonance gas model. 
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Figure 2: Details of EoS data for the p4fat3 and asqtad actions. The left figure shows the range T = 
180 - 300 MeV and the right figure shows the range T = 300 - 700 MeV and fits to the p4fat3 data. 

Data in the range T = 140 — 200 MeV lie below the hadron resonance gas (HRG) model. 
There is an upward shift (larger value of (e — 3p)) on going from N t = 6 to 8, which for some of 
the points is comparable to the difference between the N t = S data and the HRG estimate. 

The discretization effects are most pronounced in the range T = 180 — 300 MeV. We find up 
to 20% decrease in the peak on going from N x = 6 to 8 for both actions and the asqtad data lie 
below p4fat3 values by up to 15%. We find that the position of the peak remains above 200 MeV, 
as already observed in the N r = 4 data with p4fat3 action f|]. 
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At higher temperatures (T > 300 MeV) the N % = 6 and 8 data are consistent. There, however, 
are two unresolved issues. The first is possible finite volume effects at T > 400 MeV. The one N T = 
8 asqtad point at T = 400 MeV on 64 3 lattices is consistent with those on 32 4 lattices, however, 
more data are needed at T > 300 MeV. Second, when fitting the data using the expected cq + 
ci/T +C4/T form, the running of the QCD coupling requires the asymptotic behavior to have a 
g 4 variation in cq. Current data are not good enough to resolve this feature. 

Overall, our preliminary results provide a reasonably consistent picture of the EoS over the full 
temperature range and are of sufficient precision to be incorporated into hydrodynamical models. 
The current uncertainty of up to 20% at 180 — 300 MeV in the trace anomaly is expected to reduce 
to roughly 5% with inclusion of data at mi/m s = 0.05 and new simulations on N x = 12 lattices. 



5. Entropy Density Across the Transition 

The energy and entropy densities are two very useful markers of the crossover temperature as 
these are the quantities that enter into the phenomenological analyses of heavy ion collisions. 

The data for energy density £, pressure p, and the entropy density s/T 3 = (e + p)/T A are 
shown in Fig. ||. We observe a rapid cross-over in the range 175 — 205 MeV and a smaller difference 
between N? = 6 and 8, and conclude that quarks and gluons become the dominant and relevant 
degrees of freedom above 250 MeV, a region that will be probed at the LHC. 
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Figure 3: Results for energy, pressure and entropy density for the p4fat3 and asqtad actions with me/m s = 
0.1 over T = 140 — 550 MeV. The band T = 185 — 195, drawn to guide the eye, covers the inflection point. 



6. Deconfinement and Chiral Transitions 

This section summarizes results for quantities used to probe the deconfinement (Polyakov loop 
and quark number susceptibility) and the chiral (chiral condensate and its susceptibility) transitions. 

The data for the renormalized Polyakov loop (L ien (T)) = Z(g 2 ) Nr (Lbare(T)), which measures 
the free energy of an isolated quark, L ren = exp[— F 00 (r)/(7 T )], is shown in Fig. ||. There is a 
small, < 10%, difference between p4fat3 and asqtad data above T = 200 MeV with the general 
trend that the difference decreases with increasing N t . The continued slow rise beyond T ~ 250 
MeV results in a broad shoulder past the peak in the associated susceptibility, and by N t = 8 there 
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no longer is a well denned inflection point (in L ren ) nor a peak in its susceptibility. We find that 
locating the transition T from the inflection point in L ren is already marginal at N T = 6 
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Figure 4: Deconfinement indicators for p4fat3 and asqtad actions. (A) Renormalized Polyakov loop. (B) 
Strange quark number susceptibility Xs/T 2 - The band 185 — 195 MeV captures the inflection point. 

P4fat3 and asqtad data for the strange quark number susceptibility, VTXs = dh\Z/d 2 {pL s /T), 
are shown in Fig. |J They are consistent for both N t = 6 and 8 and show a rapid crossover with 
the inflection point covered by the band at 185 — 195 MeV. % s is a good probe of deconfinement as 
it measures fluctuations in the strange charge and therefore it does not need renormalization. The 
location of the peak in the fluctuations, however, requires calculating the fourth derivative, which 
is still in progress. 

The chiral condensate is investigated using the combination 

A f T ) - (W)e(T)-m e /m s {\jhir) s (T) 
l ' s (xj/xjf)i(0)-me/m s (xi/xjr) s (0) 

in which the additive ultraviolet divergence of the form m/(a 2 ) at nonzero quark mass is removed. 
The data in Fig. || show consistency between the two actions by N T = 8 indicating that the residual 
different additive and multiplicative renormalization factors due to slight differences in the lattice 
parameters between the two actions mostly cancel in A/ , S (T). The difference between N x = 6 and 
8 data are mostly accounted for by shifting the T of N t = 6 data to the left by 5 MeV. Lastly, the 
inflection point in the crossover is again captured by the band at 185 — 195 MeV. 

Fluctuations in the light quark condensate are given by the isosinglet chiral susceptibility. It 
consists of the connected and disconnected parts Xsingiet = Xdhc +2# C onn- Data for ^d isc for the 
two actions on N r = 8 lattices are compared in Fig. ||] and exhibit four features. (A) Consistency 
between the two actions; (B) an increase in peak height with decreasing quark mass; (C) A shift 
by ~ 10 MeV towards smaller T between me/m s = 0.2 and 0.1; and (D) a broadening of the 
peak as m q — > 0. Karsch |7|] has argued that, below the transition in the broken chiral symmetry 
phase, fluctuations caused by the vanishing Goldstone pion mass should broaden the peak, i.e., 
Xsingiet —> 00 with m q — > for all T < T c . If this picture (consistent with our data) is substantiated 
then the location of the right edge of this broad peak would be the appropriate locator of the chiral 
symmetry restoration temperature. Extrapolating this edge, using our m(/m s = 0.1 and 0.2 data, to 
the chiral limit gives Transition > 180 MeV for Af T = 8. 
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Figure 5: P4fat3 and asqtad data for the chiral symmetry restoration transition: (Left) The subtracted con- 
densate A/ s ; (Center) Xdisc f° r m t = 0.1m. s ; (Right) Comparison of peak in Xdisc with m f/ m s = 0.1 and 0.2. 
Note the broadening of the peak and strong m qmr k dependence attributed to Goldstone modes [j7j] . 

7. Conclusions 

Data for (e — 3p) /T 4 show consistency between asqtad and p4fat3 actions and exhibit < 20% 
change between N t = 6 and 8. We consider estimates of e and p on N t = 8 lattices are precise 
enough to use in phenomenological analyses of the evolution of the QGP at RHIC (T < 300 MeV). 

All observables used to probe the transition (energy and entropy density, Polyakov loop, quark 
number susceptibility, chiral condensate and its susceptibility) show a rapid crossover that takes 
place between 175 <T< 205 MeV with 185 — 195 MeV as our best estimate for the transition 
temperature most relevant for phenomenological studies. 

These calculations are being extended in two ways. First, simulations are being done at 
mi = 0.05m v to directly probe the system at approximately the physical u,d mass \m. Second, 
for continuum extrapolation we plan to simulate N x = 12 lattices in addition to N z = 6 and 8. 
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